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We explore the potential of direct spectroscopy of ultra-narrow optical transitions of atoms localized 
in an optical cavity. In contrast to stabilization against a reference cavity, which is the approach 
currently used for the most highly stabilized lasers, stabilization against an atomic transition does 
not suffer from Brownian thermal noise. Spectroscopy of ultra-narrow optical transitions in a cavity 
operates in a very highly saturated regime in which non-linear effects such as bistability play an 
important role. From the universal behavior of the Jaynes-Cummings model with dissipation, we 
derive the fundamental limits for laser stabilization using direct spectroscopy of ultra-narrow atomic 
lines. We find that with current lattice clock experiments, laser linewidths of about 1 mHz can be 
achieved in principle, and the ultimate limitations of this technique are at the 1 \x Hz level. 



Ultra-stable lasers are central components of optical 
atomic clocks and precision spectroscopy. Today's most 
stable lasers are made by locking the frequency of a 
prestabilizcd laser to a resonance of a high finesse ref- 
erence cavity [IH3] . The phase stability of these lasers is 
limited by thermal noise in the mirrors of the reference 
cavity @|. They achieve linewidths below 1 Hz [5] corre- 
sponding to oscillator quality factors (Q-factors) of order 
10 15 . Improving laser stability beyond the current state 
of the art will have a significant impact on precision sci- 
ence and quantum metrology [5] , but further advances in 
laser stability through refinement of reference cavities re- 
quires a significant investment in resources given the ma- 
turity of the optical designs involved [7] . The purpose of 
this paper is to propose an alternative laser stabilization 
technique, by means of direct cavity-enhanced nonlinear 
spectroscopy, and to elucidate the rich phenomenology 
of this approach in an extreme regime of cavity quantum 
electrodynamics and optical bistability. 

Strong optical transitions typically used for laser stabi- 
lization are not suitable for ultimate laser stability since 
the atomic transition frequency is very sensitive to stray 
fields, collisions, etc. However, for special ultra-narrow 
optical clock transitions that are now being routinely 
used for optical atomic clocks [8HTU]. these shifts are 
small, very well characterized, and can in some cases be 
eliminated or controlled [XT] . 

Compared to the use of strong transitions, the physics 
of this frequency locking scheme is non-trivial because 
the atomic transition is strongly saturated for very small 
intensities. Additionally, sufficient free-space optical 
depths are not available in current-generation experi- 
ments. One can circumvent this problem by working 
in a cavity-enhanced, highly non-linear, strongly satu- 
rated regime in order to achieve a signal that is strong 
enough for laser feedback. This regime has been studied 
extensively in the context of non-linear optics with alkali 
atoms [12U15] . albeit in a much less extreme limit. 

In this paper we consider a simplified model that con- 



tains all the essential components of this many-atom cav- 
ity QED system (Fig. [I]), but in the extreme bad cav- 
ity limit. Here, despite the unavoidable strong satura- 
tion effect, we are able to uncover a collective atomic 
interaction regime where we preserve the superior fre- 
quency discrimination capability of a narrow atomic tran- 
sition. This model serves as a basis upon which to cal- 
culate the fundamental limitations of our stabilization 
scheme, although real- world implementations will require 
more complicated topologies. One such approach could 
be based on the NICE-OHMS technique [M HZ], where 
the local oscillator and signal beams are co-propagated 
through the cavity to reject common-mode frequency 
noise. The effects of finite vacuum lifetime and heat- 
ing could be addressed by operating two systems in a 
multiplexed fashion, while heating could additionally be 
mitigated at the single-system level by implementing a 
Raman cooling scheme similar to that proposed in [T5] , 

In our simplified theoretical analysis, we consider an 
ensemble of N two-level atoms with transition frequency 
w Q trapped in an optical lattice potential inside a cavity. 
The lattice is at the magic wavelength where the differ- 
ence of the AC Stark shifts of both levels vanish [19] . The 
atoms are assumed to be in the vibrational ground state 
along the lattice direction and in the Lamb-Dicke regime 
such that we can neglect Doppler broadening and recoil 
effects. The atomic transition is near resonant with a 
cavity resonance with frequency lu c and field decay rate 
k. A laser with frequency ujj, is coupled into the cavity 
and the transmitted light is detected by means of bal- 
anced homodyne detection. 

This A-atom system is described by the Hamiltonian 
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FIG. 1. (Color online) Schematic of cavity enhanced ultra- 
narrow linewidth absorption spectroscopy for laser stabiliza- 
tion, M: mirror; BS: beam splitter; LO local oscillator; D: 
photodiode. 



Here, a z — \ e j)( e j \ ™ \9j){9j\ ^ s the operator for the in- 
version of atom j, and d+ — \ej)(gj\ and a_ — \gj){ej\ 
are spin raising and lowering operators, respectively. 
The bosonic field operator a is the annihilation opera- 
tor for a photon in the cavity. The coupling constant 
g — (p/h) y/Huj c / (2T4ff£o) is half the vacuum Rabi fre- 
quency with V e g the effective mode volume of the cavity, 
p the dipole moment of the atomic transition, and eo 
the vacuum permittivity. The cavity is classically driven 
with amplitude rj by the in-coupled laser. 

In addition to the coherent dynamics described by the 
Hamiltonian we also need to account for dissipative pro- 
cesses. These are spontaneous emission from the excited 
atomic state (decay rate 7), decay of the atomic dipole 
with rate T 2 _1 , and decay of the cavity field with rate 
k. We treat these dissipative processes within the usual 
Born-Markov master equation [2D]. Although we do not 
consider inhomogeneous atom-cavity coupling, this effect 
does not change our results qualitatively and can in prin- 
ciple be taken into account primarily by a rescaling of the 
cooperativity parameter via an effective atom number. 

We assume that the cavity is locked to the probe laser, 
i.e., u>l — u c . This could be achieved for example by us- 
ing a frequency-offset Pound-Drever-Hall locking scheme 
|21j on a different cavity longitudinal mode in conjunc- 
tion with a piezo-tuneable cavity. Effects due to a slight 
detuning between laser and cavity are negligible owing to 
the comparatively large cavity linewidth and we further 
quantify this statement in Appendix |B| 

To study the non-linear dynamics of this system we 
consider a semi-classical approximation where all expec- 
tation values of more than one operator can be factorized, 
e.g. (a^a^) ss (al){a_). Consequently, we find the set 
of first order equations of motion for the expectation val- 
ues o = (6) with o € {a , <r_ , a z }, 
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FIG. 2. (Color online) Intra-cavity intensity as a function 
of incoupled intensity for C = 100 and A = (blue solid 
line), A = 10T 2 _1 (purple dashed line), and A = lOOT-f 1 
(yellow dotted line). The vertical dashed lines mark the 
lower and upper threshold for bistability. The diagonal 
dashed lines show the intra-cavity intensity for a completely 
bleached atomic ensemble, \a\ 2 — Ii n , and for the unsatu- 
rated limit, \a\ 2 = 7i n /(2C) 2 . The saturation photon number 
isn = 7 r 2 "7(45 2 ). 



The atom-cavity detuning is A = ui a — lj c . 

The steady state of the system is obtained by setting 
the time derivatives to zero. The steady state polariza- 
tion of the atoms is given by 
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(5) 



Inserting this into the equations for the inversion we find 
the saturated inversion 
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where no = 7T 2 1 /(4g 2 ) is the saturation photon number. 
The mean number of photons in the cavity is then 
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Here, C = NCq is the cooperativity parameter and Cq — 
g 2 /(kTo -1 ) is the single atom cooperativity parameter. 

In this proposal we consider a regime of high- 
cooperativity where the total optical depth of the atom- 
cavity ensemble is greater than unity in the weak-driving 
limit. Specifically, in order to enter the nonlinear regime 
of spectroscopy considered here, the total cooperativity 
must satisfy C > 8. The solution for the steady state 
intensity with C = 100 is illustrated in Fig. [2] For low in- 
coupled intensity, I ln = rj 2 / (no^ 2 ) < 4C, the atoms and 
cavity behave like two coupled harmonic oscillators. For 
uj a = lu c the resonances of the coupled system are split 
by 2gi/N, the vacuum Rabi splitting. Hence, the driving 
field is far detuned from the coupled-system resonances 
for A = and the intensity inside the cavity is reduced 
by a factor 1/C 2 compared to an empty cavity. On the 
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FIG. 3. Intra-cavity intensity as a function of detuning of the 
driving laser from the atomic resonance and from the cavity 
for C = 100 and 7 in = 5 x 10 3 . Only the solution with the 
largest intra-cavity intensity is shown. Near resonance there 
are two additional solutions (see Fig. HJ. The white hyperbo- 
las indicate the resonances of the weakly driven system. 
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FIG. 4. (Color online) Intra-cavity intensity for C — 100 as 
a function of detuning for incoupled intensities Ii n = 5 x 10 3 
(blue solid line), I- m = lx 10 3 (purple dashed line), and /i n = 
1 x 10 2 (yellow dotted line). 



other hand, in the strong driving limit, /;„ > C 2 /4, the 
atomic transition is completely saturated and the cavity 
behaves as if it were empty. In the intermediate regime, 
AC < Ii n < C 2 /4, two stable solutions exist; a low inten- 
sity branch on which the atomic transition is unsaturated 
and a high intensity branch on which the atomic transi- 
tion is saturated. 

To clarify the connection of the physics considered here 
with previous studies of optical bistability in cavity QED, 
it is useful to consider the intensity in the cavity as a 
function of A and lo c — ljl ■ One of the stable solutions 
for the intra-cavity intensity is shown in Fig. [3] In the 
weak driving limit, I ln 0, the resonances of the system 
approach the white hyperbolas while the resonance of the 
strongly driven system, I ln — > oo, lies on the black hori- 
zontal line. Remarkably, with the axis rescaled as in that 
figure, the plots depend only on two free parameters, C 
and 2; n . Most experiments on optical bistability in cav- 
ity QED to date have been carried out in a regime where 
C/T2 3> k. For such an experiment, scanning Ul with 
Lj a = co c corresponds to the nearly vertical dotted line in 
this figure [23] ■ In our proposal A is scanned while 
lj c = ujl at all times, corresponding to the black horizon- 
tal line. While the basic physics behind this non-linear 
coupled system has been known for a long time |22j . it 
has not been interrogated in the way discussed here. 

The spectra resulting from scanning A in this way are 
shown for weak, intermediate (i.e., bistable), and strong 
pumping in Fig. |1J These spectra are cuts through 
the plot in Fig. [3] along the A = line. In the weak 
pumping regime (dotted line) we see a broadened ab- 
sorption feature with width CT 2 _1 . In the bistable regime 
(dashed line) there are three possible stationary values of 
the intra-cavity intensity near resonance. The solutions 
corresponding to largest and smallest intensity are dy- 
namically stable while the intermediate intensity solution 
is dynamically unstable. In the strong pumping regime 
(solid line) there is only one steady state for any detun- 
ing and a peak develops near resonance. Physically, this 



peak emerges because near resonance the atomic transi- 
tion is strongly saturated, whereas away from resonance 
the cavity field experiences an additional phase shift due 
to the atoms and does not build up in the cavity. 

In this work, the new idea is to lock the probe laser 
and cavity to this strongly saturated resonance feature. 
To estimate the potential performance of such a lock, we 
need to know the signal power and the slope of the phase 
across the resonance. The signal power is equal to the 
power leaking out of the cavity in steady-state, and is 
given by 

P ~ hu L KC 2 n f3/2 = 2huj L ri 2 /K , (8) 

where the parameter j3 = 4I nl /(C 2 ) > 1 describes how 
far above the upper threshold for bistability the sys- 
tem is driven. This power corresponds to a photon shot 
noise limited bandwidth-normalized signal to noise ratio 
of SNR 2 = nC 2 n j3 [Hz], assuming unity photodetector 
quantum efficiency. Specifically, near resonance, we can 
write the differential photocurrent from the system as 

Miff = ^ V p sigPw^5v (t) + Si (t) . (9) 

Here, <p is t nc frequency-dependent phase shift imparted 
by the intracavity atomic medium near atomic resonance, 
5i (t) is the shot noise noise on the photodetector dif- 
ference signal, 8v (t) is the system detuning from exact 
atomic resonance, and ^LO(sig) is the optical power in 
the LO (signal) pathway. Shot noise will contaminate 
the resonance condition as 

The phase shift near atomic resonance is linear to first 
order for small frequency deviations, and is given by 
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TABLE I. Quantum limited linewidth according to Eq. (A19l for several optical lattice clock systems. The cavity geometry 
is 14ff = L x (100 A™) 2 and the finesse, T, is tuned to give NCo — 100. The signal to noise ratio (SNR) is 1 Hz bandwidth- 
normalized. In all but the last case, Ti values have been set to be <1 s. This is a conservative estimate based on current- 
generation lattice clock experiments |23| . 



The shot-noise limited photocurrent noise has a white 
power spectrum and in the limit of Plo 3> -Psig, the mag- 
nitude is proportional to ^Plo- As a consequence, the 
frequency noise power spectral density of the lock error, 
Sg v , is white. We convert this quantity to conventional 
laser linewidth (see e.g., [211 [25] and Appendix |A|) when 
the system is locked and find that 
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This is the key result of this paper, as it represents the 
quantum-limited linewidth, Av, of a laser stabilized to 
the nonlinear resonance feature discussed in this work. 

It is worth contrasting these results with the ones ob- 
tained for a proposed active laser based on ultra-narrow 
optical transitions [TH]. For that system the linewidth 
is given by Az/i asor = Coj/ir. The atoms behave more 
collectively in the case of the laser. At the peak of laser 
emission the collective dipole of the atoms is proportional 
to N, i.e., (J+J-) oc N 2 , where J_ = ,P + = Ylj=i ^ n 
contrast, for the passive spectroscopy considered here 
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on resonance, A = 0, i.e. the effective number of atoms 
that participate in the collective dynamics is reduced by 
a factor of order y/T^y. Finally, we note that in the limit 
where there is no inhomogeneous broadening (T 2 — 2/7), 
Eqn. §Alty reduces to Av = f3C j/ (64tt). This is, for p 
of order unity, the same scaling as in the laser case. 

Table [I] summarizes the stabilization performance that 
can be achieved for several atomic species and transi- 
tions. In all these examples the parameters are cho- 
sen such that C « 100. The mode volume of the 
cavity is = L x 7r(100/im) 2 , where the length L 
does not enter the results. Furthermore, in this lock- 
ing scheme, the quantum-limited lock bandwidth (be- 
yond which the signal to noise drops below unity) is given 
by BW q i = nC 2 n p. In all cases considered, this funda- 
mental limitation is well above the kHz range, indicat- 
ing that the requisite level of laser pre-stabilization is 
well within current technological capabilities. In several 
realistic lattice clock systems, we find that laser stabi- 
lization can achieve quantum-limited performance at the 



mHz level without suffering from thermal noise. Finally, 
improvements in the coherence time T 2 of the narrowest 
transitions yields reciprocal gains in the quantum-limited 
locked-laser linewidth, underscoring the importance of 
investigating possible decoherence mechanisms for neu- 
tral atom lattice clocks beyond the 1 s time-scale. 

In conclusion, we have proposed a laser stabilization 
technique based on strongly saturated spectroscopy of 
narrow optical transitions that enables linewidths in the 
1 mHz range with current experimental technology. This 
technique is not limited by thermal noise and the fun- 
damental limits of this scheme are below the 1/uHz level. 
In the future we plan to study alternative realizations of 
this idea including atomic beams and trapped ions. 

We thank J. K. Thompson and J. Cooper for valuable 
discussions. This work has been supported in part by 
NIST, NSF, DARPA, and ARO. 



Appendix A: Derivation of locked laser linewidthd 

In this appendix we derive in detail the expression for 
the quantum noise-limited linewidth, which is presented 
in Eqn. 12. We begin by considereing the configuration 
shown in Fig. [T] The photocurrents of detectors one and 
two are given by 
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with the "+" ("— ") corresponding to detector one (two). 
Here, Aip is the additional phase shift acquired by the 
signal beam, Plo is the power in the local oscillator path- 
way, P s jg is the power in the signal pathway, ?7 qo is the 
detector quantum efficiency, and ^1,(2) is the stochasti- 
cally fluctuating component of the photocurrent at detec- 
tor one (two) due to shot noise. Thus, with the proper 
choice of LO phase and assuming Aip 1, 



Miff (t) = h - 12 



qc 



yfP~P^Atp + 6%i (t) - 5i 2 (t) . 

(A2) 

We re-write the term 8i\ (t) — 5i 2 (t) as Si (t) = 8i\ (t) — 
6i 2 (t). The time-domain autocorrelation of 5i (t) is given 
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FIG. 5. Phase shift of the transmitted cavity light (with re- 
spect to the LO) due to the atomic medium inside the cavity 
and as a function of laser detuning. Here, C = 100 and /3 = 2. 
The dotted line is the linear approximation for the phase near 
zero detuning. Inset: Cavity transmission curve for the same 
parameters. 



by 



(Si (t) Si (t + t)) 



e 2 ?7qe 

hv 

e 2 r) qc 
hv 



[Psi S + PLo}S(r) 
PwS (t) . 



(A3) 



Here, S (r) is the Dirac delta function. This corresponds 
to a two-sided photocurrent noise power spectral density 
of 
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The resonance center is observed in this system via 
the difference photocurrent — namely where the difference 
photocurrent is equal to zero. In order to see the effect of 
the LO shot noise on the lock stability, one can expand 
idiff about zero detuning to linear order of Av as 



Miff = ^VPsigP L oA</? + Si (t) , 

c 2 -^^rp^^Av + Si(t). 

hv av 



The validity of making this linear approximation is shown 
in Fig. [5j where the complete phase shift of the medium 
is shown as a function of detuning from resonance along 
with an analytical solution for the linear phase shift, 
^Av. As long as the laser is close to resonance, the 
phase is linear to a good approximation. Now we can 
identify the effect of Si (t) on our ability to determine 
the line center of the atomic resonance. 

Under locked condition, the DC difference current, 
/diff, is enforced to be zero via control of the laser fre- 
quency. We can thus see that the term Si (t) corrupts 
the measurement. Namely, our frequency error is given 
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The denominator came directly from Eqn. |A5| 

When locked to the cavity-atom resonance, we assume 
the laser has an electric field given by 

E(t) = £ oe i2™ t+Wt). ( A7 ) 
Here the phase error, S<f)(t), is related to Av crr (t) by 

^ = Av CII (t) . (A8) 

In order to go from this time-domain expression to the 
frequency domain via the Wiener-Khinchin theorem (fol- 
lowing the general approach presented in Ch. 3 of |26|). 
we compute the autocorrelation of the field amplitude, 
Re (t), given by 

R E (t) = (E(t)E*(t + T)) 

= |£ |V 2 ™ or (e lW(t) ~^ ( * +T)) ). (A9) 

An application of the gaussian moment theorem gives 

( e i(m)-mt+r))} = exp {-([^(t) - 6<f>(t + r)] 2 ) /2}. 

(A10) 

We can re-write the expectation value as 

([ty(t) - 5<p{t + r)f) = 2([^(r)] 2 ) - 2([S<f>(t)S<f>(t + r)]> 
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It is then a direct consequence of the Wiener-Khinchin 
theorem that 
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where Ss,j, (/) is the two-sided phase fluctuation power 
spectral density for 8<p(t). However, we can easily re- 
late Ss<p to S'ai/ (the two-sided frequency deviation power 
spectral density) by Eqn. A8 such that 
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Applying the Wiener-Khinchin theorem to Eqn. |A6| 
we have 
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We can therefore re- write Eqn. |A13| as 
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with ho given by 



input/output dynamics of the system 
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Now we have an expression for the electric field autocor- 
relation, namely 

R E ( T ) = \E \ 2 ( Al7 ) 

We apply the Wiener-Khinchin theorem to this expres- 
sion and obtain a Lorentzian profile for the laser op- 
tical power, with frequency full width half maximum, 
Ai^FWHM, given by 
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We combine this with the results of Eqns. 8 and 11 of the 
main text, and obtain the result presented in Eqn. 12, in 
the limit of unity detector quantum efficiency, namely 
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y = x\ 1+C- 
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The parameter x is related to (a) by x = (a) / ' ^/n^, y is 
given by r\j (fty 7 ^)) (|y| 2 = I[ n ), A and T 2 are the same 
as given in the text, and the parameter 9 is the cavity- 
laser detuning in units of n, namely 9 — (co c — lji) / k. In 
the text, 9 was assumed to be negligibly small. Here we 
quantify this statement. 

If we assume that we are near resonance in the non- 
linear, strongly saturated regime (/3 > 1, C > 1), then 
\y\ 2 ~ \x\ 2 = l3C 2 /4. If T 2 A < C, then we can expand 
Eqn. [BT] such that 



y ~ x 
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Therefore, the phase shift of the transmitted light is given 
by 



A</> = Arg [x/y] 
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From Eqn. |B3| it can be seen that for a given cavity- 
laser detuning, the lock center frequency shift, A^i ascr , is 
given by 



A leaser — 



C/3 fui c - lji 
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Cavity lock precisions of > 10 4 are routinely achieved in 
the laboratory. This implies that 



Appendix B: Line-pulling effects due to cavity-laser 
detuning 



In order to derive the line-pulling due to an imperfect 
lock between the cavity and probe laser, we make use 
of the full optical bistability equation that describes the 



< 10" 



(B5) 



For typical parameters considered in the main text, 
namely C = 100, p = 2, and T 2 = 1 s, this implies that 
the cavity pulling effect is below the 1 mHz level. Longer 
T 2 times will further suppress this effect. 
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